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A limit on the tau neutrino mass is obtained usins all the Z° t+t' data collected at LEP by the DELPHI 
detector between 1992 and 1995. In this analysis events in which one of the taus decays into one charged particle, 
while the second t decays into five charged pions (1-5 topology) have been used. 

The neutrino mass is determined from a bidimensional fit on the invariant mass ml^ and on the energy E^-^ of 
the five TT* system. The result found is m^^ < 48.0 MeV/c^ at 95% confidence level. 



1. Introduction 

The question whether neutrinos are massive is 
one of the outstanding issues in particle physics, 
astrophysics and cosmology. 

Among the possible frameworks, the "see-saw" 
mechanism is considered to be the most inter- 
esting because it explains the smallness of neu- 
trino masses by connecting them to the scale of 
new physics. In fact the see-saw mechanism as- 
sumes a neutrino mass hierarchy similar to that 
of quarks or leptons. 

On the basis of cosmological arguments a sta- 
ble tau neutrino with a mass larger than a few 
eV/c^ cannot exist; however unstable neutrinos 
can be more massive, so a direct neutrino mass 
measurement of the order of a few MeV/c^ would 
suggest that the Vr is unstable 

The clearest evidence of neutrinos having mass 
comes from the oscillation experiments ||^ , even if 
this kind of experiments cannot give the neutrino 
mass scale. Only direct measurements can do it 
and the most recent experiments [Q have fixed 
the mass of the Ve to be about 2 eV/c^. However 
the results of the oscillation experiments must be 
confirmed before concluding that massive unsta- ^ 
ble neutrinos cannot exist. 

The kinematic properties of the observed final 
state in hadronic r decays give direct bounds 
on TOj^^ . The best decay channels for a missing 



mass measurement are those with a high hadronic 
mass, where the little energy is available for unde- 
tected neutrinos. This kind of r decays presents 
a small branching fraction , so the data sample 
generally consists of a few events. 

High multiplicity decays 3n^2n^T7r{i'T) 
have been selected for this analysis. The end 
point of the invariant mass distribution of the 
hadronic system, compared to rur, gives a limit 
on the Vr mass. 

Several direct experimental measurements have 
been done on the Vt mass. The world's best limit 
was obtained by the ALEPH Collaboration 
with nit,^ < 18.2 MeV/c^ at 95% confidence level. 

2. The method 

Multi-hadron decays of the r lepton |^,|| are 
due to the couplings of the charged weak current 
to hadrons, which in the case of bn^Vr 
decays reduces to the coupling of the boson 
to the ud' current. 

We can describe the tau decay as a two-body 
decay: 

{Er,Pr) -> h- {Eh,ph) + Vr {E^,Pt.) (1) 

where the hadronic system h~ is composed by 
five pions. In the tau rest frame, the energy of 
the hadronic system is used to compute the value 
of m^^ because of its dependence on nir , ruh and 
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In the laboratory frame this energy becomes: 



(2) 



Eh ^ {El + PpI cost 



(3) 



where (5 is the tau velocity, 7 is the boost factor 
and is the angle between the direction of the r 
and that of the hadronic system in the tau rest 
frame. We assume that the r energy, E^, is equal 
to the beam energy, Eiieam ■ Initial (ISR) and final 
(FSR) state radiation can reduce Er and will be 
considered later in section |^. 

Because of the undetected neutrino, the direc- 
tion of the T is not determined, so rriij^ cannot 
be computed directly. However, the energy of the 
hadronic system is helpful because it partially re- 
covers the loss of information. In fact Eh de- 
pends on cos 6 and it must fall inside the interval 

Emin — Eh ^ E^f-^^x • 



E 



min.max 



^{El±(3pl) 



(4) 



As a consequence kinematic allowed regions are 
defined for different values assumed by m^^ , as 
shown in Fig. 0. 

The probability to find an event in the allowed 
region clearly depends on the kinematic of the 
decay and on the resonance structure. Because 
both the invariant mass mg^ and the E^^^ / Eheam 
of the hadronic system are functions of the mea- 
sured momenta, a positive correlation arises be- 
tween these two quantities. The error ellipses for 
each event must be taken into account, so that, 
as shown in Fig. ^, event 1 constrains the neu- 
trino mass m^^ much more than event 2, even if 
the latter has a higher hadronic mass. This is the 
advantage of fitting the distribution of the two 
variables and i?57r/i?6eam, rather than 
alone. 

3. The DELPHI detector 

The DELPHI detector and its performances are 
described in detail in [p|jic|]. The whole detec- 
tor, composed by a main central part (barrel) and 
by two forward regions (end-caps), covers the full 
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Figure 1. Two hypothetical — > Btt^Vt events 
have been plotted in the (ml^ , E^^^ / Ebeam) plane 
with their error ellipses. The coloured lines bound 
the kinematic allowed regions for different values 
of rrii, . 



solid angle. A superconducting solenoid provides 
a uniform axial magnetic field of 1.23 T. Charged 
particles are tracked in the barrel region using a 
combination of four cylindrical subdetectors: the 
silicon Vertex Detector (VD), the Inner Detector 
(ID), the Time Projection Chamber (TPC) and 
the Outer Detector (OD). Electron and photon 
idenification is provided by the High density Pro- 
jection Chamber (HPC) located outside the OD. 

4. Data selection 

For LEP centre-of-mass energies of ^/s ~ M^o , 
and T~ are produced back-to-back (ignoring 
radiative effects). Each r decays producing one, 
three or more charged particles in addition to one 
or two neutrinos and, possibly, neutral mesons. 
All particles apart from neutrinos can be detected 
by DELPHI, e+e" Z" ^ t+t' decays are 
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easy to identify because the two back-to-back 
jets are narrow and with low charged multiplicity. 
The presence of neutrinos implies that not all the 
energy in the event is seen and the invariant mass 
of the final state particles is less than TOt- 

The background from multi-hadron production 
e+e" — > qq a.t LEP energies can easily 

be reduced since it is characterized by a rela- 
tively high charged multiplicity and by a high in- 
variant mass. Final states from e^e~ — > e^e~ 
and e+e^ — > A^^A*^ leave very characteristic sig- 
natures in DELPHI and are effectively discrimi- 
nated. 

The analysis presented here is based on the 
data collected by DELPHI from 1992 to 1995; the 
centre-of-mass energy y/s of the e"'"e~ system was 
between 89 and 93 GeV. 

4.1. e+e~ — >■ Z° — !• x+T~ selection 

The first step in data selection was to reject 
all non— Z'^ decay (cosmics, secondary interac- 
tions in the vacuum pipe, etc.). Each event has 
been divided into two hemispheres defined by the 
plane perpendicular to the thrust axis. It was re- 
quired that the reconstructed topology was one 
charged particle in one side and five in the other 
(1-5 topology). 

Events were accepted if Pms, the sum of the 
momenta of all the charged particles plus the to- 
tal neutral electromagnetic energy, was greater 
than 0.09-v's, rejecting some of the hadronic 
decays, 77 interactions and beam gas interac- 
tions. To reject cosmics and to reduce further 
the contamination from beam gas events it was 
required that, defining particles 1 and 2 as the 
most energetic in each hemisphere, ri^2 < 1-5 cm 
and \zi^2\ < 4.5 cm, where rj and zj are their 
transverse and longitudinal distances of closest 
approach to the average beam spot. 

Part of the hadronic qq background was re- 
jected by the requirement on the topology and on 
the impact parameters as described above. The 
remaining low multiplicity qq and 77 events can 
be easily rejected by requiring that the isolation 
angle 9iso be greater than 160°. The isolation an- 
gle is defined as the minimum angle between any 
pair of tracks belonging to opposite hemispheres. 



Table 1 



Composition of the selected data sample 



Background Class 


Back^ 


ground Source 




5^±^o 


8.3% 




3^± 


2.3% 




37r±n7r° 


3.4% 


non— T 


qq 


all at ml^^ > rur 



(*) in the fit region it is reasonably flat 



Another discriminating variable was the total 
momentum of the 5-prong system, P5. It must be 
greater than 30 GeV/c. Even if either 0iso and P5 
are related to the invariant mass of the hadronic 
system, the distribution of 7715^ in the region of 
the end point is not affected by these selections. 

To reduce further qq it was required that the to- 
tal number of neutrals must be less than six. The 
hadronic decays usually have a large number 
of neutrals, whereas r+r^ events should have no 
or few neutrals. Four-fermion background was 
suppressed by the requirement on the topology. 

4.2. 37r=t 277^17.^ (/y^) selection 

The minimum momentum among the five 
charged particles and the minimum number of 
hits in the Vertex Detector help to reject those 
events that are reconstructed with a 5-prong mul- 
tiplicity because of a secondary interaction or a 
tt'^ Dalitz decay (tt'^ — > 7e+e~), as for example 
Sn^i^r > 7r°. For this kind of events, sec- 
ondary interactions or successive decays cause a 
degradation of the outcoming momenta. Events 
with Pmin less than 0.5 GeV/c were rejected. 

The minimum invariant mass of pairs of op- 
positely charged particles in the hemisphere con- 
taining five tracks, m*g (assuming the electron 
mass for both), was required to be greater than 
0.05 GeV/c^, to reject events with a 7 ^ e+e~ 
conversion. The bn^i^r > ""^ decays are 

rejected requiring the maximum value of the elec- 
tromagnetic energy being less than 4 GeV. 

4.3. Final data sample 

The composition of the selected data sample, 
as determined with the Montecarlo, is listed in 
Table |l|. 
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Figure 2. Distribution of the hadronic invariant mass from simulation {left) and the comparison between 
Montecarlo and data {right). 



A total of 47 5n^i/r candidates have been 

selected in the data. Fig. || shows the distribution 
of the hadronic invariant mass mt^^ for the simu- 
lation (signal and backgrounds) and a comparison 
between the simulation and data. qq events 

have values of the invariant mass greater than 2 
GeV/c^; other kind of r decays have small values 
of rUg^, so they are located on the left side of the 
invariant mass spectrum. 

5. The likelihood function 

In order to determine the upper limit on m^^ an 
unbinned likelihood fit has been performed. The 
likelihood function depends both on the invariant 
mass and on the energy E^t^ of the hadronic 
system: 



£=11 P{m*,E,\m,^) 



observing each selected event i at the position 
(to* = TOg^, Ei EE E^^/Ebeam) in the kinematic 
plane defined as a function of the tau neutrino 
mass TOj/, : 



P{mlE,\m,^) 



(5) 



where A^obs is the number of the selected candi- 
dates, while P {m* , Ei\m^^) is the probability for 



1 d^N{mlE,\m,^) 
' N dm dE 
'TZ{m,E,p,am,crE) 

I £ (to, E) 



(6) 



where TZ is the experimental resolution function 
of the detector and e is the selection efficiency. 
Both TZ and e could be a function of either the 
invariant mass and the normalized energy. The 
shape of the resolution function TZ has been de- 
termined from simulation. Since the two variables 
are correlated, 7?. is a 2D gaussian depending on 
the three parameters am, cte and p. The invariant 
mass resolution is about 16 MeV/c^ (in average) 
and the normalized energy resolution is of the or- 
der of 10^2 ^ MeV). If the resolution 
is constant, the experimental resolution function 
can be factorized, and the likelihood C becomes 
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Figure 3. Distribution in the upper part of the 
{ml^, E5T,/Ebeam) plane for t"^ ^tt^v^ candi- 
dates. The coloured lines bound the kinematic 
allowed regions for different values of rriu^ . 



simpler. In this analysis each event is considered 
with its proper resolution. 

The selection efficiency e was found to be con- 
stant with respect to the invariant mass and the 
energy if > 1.6 GeV/c^ and E^^/Eteam 

greater than 0.85. As a consequence, the fit will 
be restricted to this region of the (m^^, i?5^/ 
Ebeam) plane. Fig. ^ shows this upper part of 
the plane with the 15 candidates used in the fit. 

The theoretical prediction P{m*,Ei\m^^) has 
been calculated as a function of m^^ including ini- 
tial (ISR) and final (FSR) state radiation with the £ 
KORALZ 4.0 ||l|] Montecarlo event generator for 
e+e" ^ Z° ^ T+T-. The KORALZ 4.0 simu- 
lator incorporated the TAUOLA 2.4 package 
for modelling the r lepton decays. 

The 1 /N factor in equation (||) ensures that the 
probability density used in the event likelihood is 
normalized for any values of neutrino mass be- 
cause the distribution of events in the kinematic 
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Figure 4. Log-likelihood of bir^Vr data 

fit as a function of the tau neutrino mass. The 
result gives m^^ < 48 MeV/c^ at 95% CL. 



allowed region of the (ml^, -Estt /-Ebeam) plane de- 
pends on the value of m^^ . 

The fraction of backgrounds has been deter- 
mined from simulation. As shown in Fig. 
hadronic events have values of the invariant mass 
greater than rriT, while background from r — > 
Stt^ > Itt^Vt or —> bTT^TT^Vr glvcs contribu- 
tion in the low invariant mass part of the spec- 
trum. The end point of rrig^ is therefore reason- 
ably background-free. However a flat background 
(6% of signal) has been considered, so the likeli- 
hood was expanded as: 

a C signal + (1 - a) Cbgd (7) 

where a is the background fraction. Only 15 of 
the 47 selected candidates are in the fit region. 
The fit gives an upper limit on the neutrino mass 
of m^^ < 48.0 MeV/c2 at 95% CL. Fig. | shows 
the log-likelihood for these 15 events. 

A possible bias in the analysis method has been 
investigated by fitting high statistics Montecarlo 
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Table 2 

Comparison between input values for the mass 
and best fit results 



m^^ (MeV/c^) 


Best Fit (MeV/c^) 


0.0 
30.0 
60.0 


7 q+30.7 
' "^-7.3 

■^^■^-32.0 

C-r 9 + 19.0 



Table 3 

95% CL m^^ limits for different decay structures 





m^^ (MeV/c2) 


CL 


pure phase space 


< 48.0 


95% 




< 47.9 


95% 




< 46.3 


95% 




< 44.5 


95% 




m * (GeV/c") 



samples with three different input values of m^^ : 
0, 30 and 60 MeV/c^. In all cases the best fit re- 
sults are in good agreement with the input values. 
The corrisponding values are listed in Table 

The expected best limit for m^^ — MeV/c^ 
is about 30 MeV/c^. 

6. Systematic effects 

Four sources of systematic errors have been 
considered. For each source a new fit has been 
performed after having changed the parameter of 
that particular source by one standard deviation. 
The difference between the 95% CL upper limit 
on the tau neutrino mass computed in section ^ 
and the one obtained by the modified likelihood 
has been considered as the systematic effect due 
to that particular source. Then all these varia- 
tions have been summed in quadrature to obtain 
the global systematic error which can be added 
linearly to the upper limit of the original fit. 

The systematic sources considered in this anal- 
ysis are summarized in the following: 

• the tau mass rur and its energy Er] 

• detector effects such as the hadronic mass 
scale and the experimental resolution func- 



Figure 5. Distributions of the hadronic invariant 
mass for data (dots) and four models of decay. 
The black hatched histogram indicates the pure 
phase space model, while the red, green and blue 
lines are the distributions obtained by means of 
three different intermediate resonance structures. 



tion. The DELPHI hadronic mass scale 
studied using the K^tt^-k^it^ decay 

shows a shift of -h4 MeV/c^ ||; 

• the selection efficiency and the effect of the 
size of the fitted region; 

• the intermediate structure of the decay. 
The hadronic invariant mass depends on 
the dynamics of the decay. Intermediate 
resonances have the effect of shifting the 
shape of mg^ spectrum to higher val- 
ues. Some studies have been performed 
by assuming a^TT^TT^^Vrivr) decay, 

— > p27r^7r^l7r (i^r) Or — > 2p7r^I7r (i^r). 

Fig. ^ shows the invariant mass distribu- 
tion compared with those predicted by a 
pure phase space model and by the mod- 
els with intermediate resonances. A model 
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Table 4 

Systematic variations of the 95% CL upper limit 
on 



Source 



Variation of m^^ 
limit (MeV/c2) 





0.3 


^beam 


<0.1 


experimental resolution 


2.5 


mass calibration 


1.0 


fit region 


0.7 


total 


2.8 


hadronic mass scale 


+4.0 



with a 2p7r^ seems to be preferred. How- 
ever, because of the small number of ob- 
served events, it is not possible to state 
which is the correct resonance. Table |^ 
shows the 95% CL upper limits obtained 
for the resonance structures considered. 

Table ^ summarises all the variations consid- 
ered in this analysis. 

7. Review of the existing collider results 

There are some interesting limits on the Vr 
mass obtained by collider experiments. Recent 
results are summarized in Table ^ and shortly 
presented in the following sections. 

7.1. Low energy measurements 

The main difference between measurements 
performed at different y/s values is that the se- 
lection between high multiplicity r decays and 
low multiplicity hadronic events from e^e" qq 
becomes easier and easier with the increasing of 
the center-of-mass energy. So LEP experiments 
select r+r^ events more efficiently and with a 
higher purity. On the other hand, taus are pro- 
duced more abundantly at CESR. 

CLEO collected tau pairs produced at a center- 
of-mass energy of 10.6 GeV. The analysis was per- 
formed using bn^Vr events and fitting the 
36 selected events, using the invariant mass and 
the hadronic energy. The 95% CL upper limit 



was m^^ < 33.9 MeV/c^ |l|. 

In addition CLEO analyzed — > 37r^27r"zyT- 
decays. This channel can be studied better at low 
energy machine, because at LEP photons from 7r° 
conversions are highly coUimated and the selec- 
tion more difficult. An upper limit at the 95% 
confidence level of 35.9 MeV/c^ was found. The 
combined result found by CLEO is 30.0 MeV/c^ 

An analysis with — )■ 37r^7r*'i/^(z/T-) decays 
was also done [|l5| . This is a useful decay because 
of the larger branching ratio than 5-prong, and in 
addition the invariant mass is very close to rur. 
From a data sample of 29058 candidates, the ob- 
tained 95% CL upper limit on the tau neutrino 
mass is 28.0 MeV/c^. 

Fig. 1^ shows the distribution of the 5 body 
selected events in the (mg^, E^j^/Ebeam) plane. 



ARGUS used r^" 



57r^ 



events collected at 



center-of-mass energy of about 10.6 GeV. Per- 
forming a fit only to the hadronic invariant mass, 
the resuh was m^^ < 31.0 MeV/c^ ||l§. 

7.2. High energy measurements 

At LEP energy, r+r^ events present a typical 
two back-to-back narrow jets structure. ALEPH 
and OPAL, as DELPHI, have analyzed data col- 
lected at center-of-mass energy equal or close to 
the Z'^ resonance. 

ALEPH obtained two separate limits by fitting 
the distribution of the visible energy and the in- 
variant mass. The fits to the 2939 — > Sn^i^r 



and to the 55 



^0 



95% CL upper limits on 



i^r events have given 
m^^ of 26.5 MeV/c2 
and 22.3 MeV/c^ respectively. These two results 
have been combined to obtain a 95% confidence 
level upper limit of 18.2 MeV/c^, and this is the 
world's best hmit on m^^ Fig. shows the 
distribution of events in the upper part of the 
(mg^, E^Tr/Eteam) plane for the 5-prong data. 
Also OPAL determined an upper limit for m^^ 



with the decay r 



57r±i/^. A limit of 43.2 



MeV/c^ at 95% CL has been obtained using a 
two-dimensional method in the hadronic invari- 
ant mass and energy distributions from 22 se- 
lected candidates |l^ . Combining this result with 
opal's previously published measurement using 
r+r^ — > Sh^V-r + ih^Ur decays, a new limit of 
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Table 5 



Limits on the 


mass in MeV/c^ from multi-hadron r 


decays 






Experiment 






Mode 












57r(7rO)i/^ 


Combined 


ARGUS 








3L0 




CLEO 




28.0 


35.9 


33.9 


30.0 


ALEPH 


26.5 






22.3 


18.2 


OPAL 


35. 3(*) 






43.2 


27.6 



(*) from T+r ih^Vr + 'ih^Vr decays 




1.66 1.68 1.70 1.72 1.74 1.76 1.78 1.80 

M.,,(GeV/c^) 



Figure 6. The hadronic scaled mass energy versus 
mass distribution of the Stt^ (a) and 37r='=27r° (6) 
data sample in the fit region from CLEO Colla- 
boration. 



m^^ < 27.6 MeV/c^ 



CL) was obtained |18| 



Before the 5-prong analysis presented here, 
DELPHI has set a limit with a measurement of 
the leptonic branching fractions of the r lepton. 
The 95% confidence level limit found was 66.0 
MeV/c^ This method is however limited by 
the precision of the leptonic branching ratio mea- 
surements, and is not competitive with the one 
described in this talk. 



8. Conclusions 

The CLEO experiment had analyzed a very 
large statistics sample and its mass and energy 
resolutions were better than LEP experiments. 

Looking at Fig. ^ and 0, it's easy to outhne 
that with comparable mass and energy resolu- 
tions, the position of events in the (m^^, E^^r/ 
Ebeam) plane determines a "/mcA; factor^'' for an 
experiment. So, the ALEPH limit is strongly 
driven by the two candidates lying in top right- 
hand corner of the kinematic allowed region. The 
most sensitive region of that plane for OPAL and 
DELPHI is empty. OPAL and DELPHI'S 5-prong 
analysis present very similar results. 
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